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Abstract 



The infrared divergence of the self-energy of a color charge is due to an enhancement 
of the long wavelength modes of the color Coulomb potential field. There are also 
long wavelength contributions to the QCD vacuum energy that are similarly enhanced. 
Vacuum modes of Hubble scale wavelengths may be affected in a cosmological setting and 
this can lead to a residual positive energy density of the form H^A^J^. Lattice studies 



constrain d. If the dark energy takes this form then the universe is driven towards de 
Sitter expansion, and we briefly study this cosmology when d is just slightly above unity. 

1 Introduction 

There have been recent suggestions [H [21 13] that in an expanding universe with Hubble pa- 
rameter H, QCD can produce a dark energy density of order HAqq^). This is an intriguing 
claim since numerically this energy density is quite close to (although somewhat above) the 
observed value, and it leads to a cosmology not that dissimilar to the cosmological constant 
cosmology. But there are reasons to be skeptical. Given that QCD has a mass gap there 
would seem to be a clear decoupling between the vastly different QCD and Hubble length 
scales. One can use the Casimir effect as an analogy, where the effects of massive fields are 
exponentially suppressed. And even with a massless field the expected contribution to vacuum 
energy is only of order H'^, again in analogy with the Casimir effect. 

For a result different from the Casimir effect, one can consider an infinite volume and simply 
remove the contribution of the large wavelength modes to the zero-point vacuum energy. Then 
the change in vacuum energy density is 
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where we give the example of a massive field with m ^ H. This could be considered the 
maximum effect that Hubble scale physics could have on a Hubble scale contribution to 
vacuum energy. But this is still a miniscule contribution and it is of the wrong sign. The 
HAqqy) behavior is said to rely on various properties of QCD, including strong interactions, 
topological effects and ghosts that decouple in flat spacetime. But it is presently unclear to 
us how the arguments of [H [2l |3] and the related work in lower dimensions [H [5] overcome the 
naive expectations. 

We shall deflne Apvac to be the H dependent contribution to the QCD vacuum energy 
density due to effects on Hubble scales. Apvac should be flnite and its dependence on H 
means that it can be discussed independently of the cosmological constant problem. From 
the previous remarks, if Ap^ac is to be large enough to be of any interest there must be some 
enhancement of the long wavelength contributions to the QCD vacuum energy. We shall argue 
that existing lattice studies can help to establish whether or not such an enhancement exists. 

We flnd particularly relevant a set of studies that are helping to elucidate the mechanism 
of confinement. The Coulomb gauge picture of confinement, originally due to Gribov [6] and 
developed over the years [^7j, connects confinement to the nontrivial structure of the gauge field 
configuration space (as manifested by Gribov copies). It has received support from numerous 
lattice studies, with [8] being a recent example. Intrinsic to this picture of confinement is the 
nontrivial scaling behavior of various quantities in the far infrared. It is this type of effect 
that is of interest to Apvac- 

On the lattice some finite volume analog of Ap^ac could in principle be measured without 
choosing a gauge. But this would be difficult in practice and so more useful is the fact that 
the lattice can provide information on the scaling behavior of 2-point functions in a particular 
gauge. This information can then be carried over to a continuum description where a gauge 
choice is necessary, and where a direct estimate of Ap^ac (a physical quantity) is more feasible]^ 
The choice of the Coulomb gauge is also appropriate for another reason, its non-covariance 
matches the existence of a preferred frame in the cosmological setting. 

What is of particular interest in Coulomb gauge is the infrared behavior of the longitu- 
dinal component of the color electric field, the color Coulomb potential. The massless long 
wavelength modes of this field are enhanced, and this enhancement must be sufficiently large 
to be consistent with confinement. We shall follow a treatment \XQl that derives a necessary 
condition for confinement in terms of the infrared divergence of the Coulomb self-energy of a 
color charge. In contrast the propagation of the transverse modes of the gauge field are sup- 
pressed in the infrared. This behavior is also related to confinement and it can be described 
by an effective mass function m{k) ^ Aqq^/A; for these modesj^ 

^The use of gauge- fixed lattice results in this way has recently been extolled in [5]. 

^It was pointed out in [5] that this would enhance the contribution in ([I]), but not sufficiently and without 
correcting the sign. 
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We then turn to the vacuum energy and explore an interplay between the nonperturbatively 
enhanced long wavelength longitudinal modes and the modes that have developed a mass 
gap. The result is an infrared enhancement of Apvac- We end section 2 by making use of a 
recent lattice study [H] to help determine the enhancement. In section 3 we return to the 
cosmological setting and discuss how Apvac feeds back and affects cosmological evolution. An 
interesting cosmology arises. 

Apvac depends on H by definition and implicit in our picture is that there is no H indepen- 
dent contribution to the vacuum energy; only then is the Minkowski vacuum consistent. Why 
this should be the case is the original cosmological constant problem. We only remark here 
on what we feel to be a necessary condition for a solution: the ultimate theory should contain 
no explicit dimensionful parameters and all masses should arise dynamically. The quantity 
Apvac of interest here is an example of a vacuum energy that appears due to the introduction 
of the dimensionful quantity H. 



2 The enhancement 

In Coulomb gauge where VjA° = 0, the color electric field of QCD can be written as -Ef — V^v?" 
where Vj-Ef = 0. Gauss's law determines the color Coulomb potential (p"" in terms of the 
transverse gluons and quarks. 

' d'yM:y{y)=p\x) (2) 

Mty = i-S^'V' + gr'^A'r{x)Vi)6{x - y) (3) 

pa ^ gq^T^q - gP^^A^EI (4) 

The Coulomb energy is 

H^oni = -I f dM^\y)W^v\y)) (5) 



2 
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d'xd'yd'z{p\x)[M-X[-^l]W-'Zp%z)) (6) 



M is the Fadeev-Popov operator and we consider its eigenfunctions in an arbitrarily large 
volume, 

j d^yMfy^iy) = X^rnix), (7) 

normalized such that 

j d'yct^ZivWAy) = Smn. (8) 

Then we have the representation 
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which gives 



mn 



where 



Fmn = I rf^xC(^)(-V^)C(^)- (li: 



This is a useful formulation since the eigenfunctions and eigenvalues of M can be directly 
determined on the lattice [101 lEl [H] • 

A diagnostic for confinement can be obtained by making the replacement 

{p'^{x)p'^{z))^^6{x-z)6ac. (12) 
This will yield the self-energy of a point color charge Q. This Coulomb self-energy is 

EseU = ^Y.^J (13) 

where F„ = In the large volume limit the sum can be replaced by an integral upon 

introducing a density of states [TO] . 



E..u = Iq'I dXpiX)^. (14) 

In the weakly interacting limit the eigenstates (p^i^) ^ire plane waves and p(A) = y/X/An'^ 
and F{X) = X, which means that itself has an ultraviolet divergence but no infrared divergence. 
But confinement requires that £^seif does exhibit an infrared divergence, and so we know that 
QCD must produce an infrared enhancement of p(A) and/or F{X). If we write 

then we see that the existence of confinement implies that a + b > 1/2 [TU] . We are using 
A = Aqcd to represent the scale below which confinement physics enters and in the following 
it will also represent the mass gap of the theory. 

Let us return to ([To| and consider the vacuum energy density pvac = -f^coui/V^- Now we 
need the quantity {p"'{x)p'^{z)) . The nonperturbative contribution to this quantity is easier 
to determine than since it is almost completely determined by the mass gap. The fields 
in p"(x) from ^ are the transverse gluons and the light quarks (of number Ug). The quarks 
develop a dynamical mass (we ignore current quark masses) and the transverse gluons have an 
effective mass as mentioned above. But the fluctuations of these fields do not in turn produce 
a mass for the Coulomb potential field. From this and dimensional analysis, the lowest order 
term in the derivative expansion of {p°'{x)p^{z)) should take the form 

{p\x)p\z))^KVl5{x-z)5ac. (16) 
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This argument does not determine the sign. 

To see more clearly where such a result comes from it is instructive to consider the one 
loop contribution evaluated in momentum space, 



k 



2 



- 1 

k 



(17) 



We have kept only the lowest order term of an expansion in p^. A now appears in ul = k'^ + A'^ 
as an approximation to the mass of both the quarks and the transverse gluons. We have 
performed a subtraction by the zero mass integral to isolate the contribution that vanishes 
with the mass gap. 

The framework here is non-covariant old fashioned perturbation theory (OFPT). The 



result in (17) is closely related to the transverse gluon and quark contribution to the usual 
vacuum polarization [13]. That contribution occurs at second order in OFPT where the 
quantity X]n(P"(^)l'^)("'l/''^('^))/(-^o ~ En) appears. The extra energy denominator turns the 



linear divergence appearing before subtraction in (17) into the standard log divergence. The 
transverse gluons and quarks give screening contributions to the vacuum polarization, and 
so the sign of our result has the same origin. Meanwhile the anti-screening contribution to 
the vacuum polarization enters at first order in OFPT and it comes from the perturbative 
expansion of 



The result of the k integral in (17) is — A7r/2, and along with the approximations = 1 
and = 3 we arrive at 

{p%x)p%z))^^AVl6{x-z)Sac. (18) 
While this perturbative derivation of a nonperturbative effect is no doubt an oversimplification. 



we shall continue to use ( 18 ) in the following. The calculation clarifies the role of both the 



mass and the screening effects of the transverse gluons and quarks. 



Inserting this result back into pvac = Hcoui/V from (10) gives 

3 A V ^ -^Ti .m. 



p ~ \ ^ ^ nm-'^ mn (19) 



mn 



According to a lattice study ^Hj the off-diagonal elements are small, so F^n ~ Fn^mn- If "we 
use this then 



Pvac ~ 



3A V — \ F?! 

-Sl/ETf (20) 



iV ^ A? 

n 

\2 



3A f^,^,,,F{X) 



dXp{X)^ (21) 



3A r .. ^/A2^('^+''') 



327r2 / V A 



dXVX — (22) 
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We are interested in the infrared behavior of this integral. Suppose there is some large 
length scale L beyond which the physics changes and the infrared enhancement effects are re- 
duced (see below). Corresponding to L would be a minimum Amin below which the integration 



in (22) would change. The modified integrand is reduced for example if it scales with a higher 



power of A for A < Amin- In any case we shall take the shift in the vacuum energy density due 



to this change to be of order the negative of the to Amin part of the integration in (22). 
To go further we need to specify how Amin scales with L, 



L2 {Lky 



(23) 



We shall interpret Amin as the minimum nonzero eigenvalue of M that exists in a finite volume 
Lr". The resulting finite shift in vacuum energy density is 

An ~ ^ A 2a+4b+l X f-°-2fe 



327r2 



mm 



« 3^A"(i)^ (24) 

where 

d={2 + c)(^-a-2h). (25) 

In the weak interaction limit a = 6 = c = 0, and d = ?>. Although it doesn't make much sense 
to consider this limit without also setting the mass gap to zero, one can see the similarity 
between ci = 3 and the naive estimate in ([T]). 

^Pvac is enhanced as d decreases, and we already see that the bound a + h > 1/2 from 
confinement will tend to do this. To do better we shall use the results of an SU (3) lattice 
study [11] 1^ There the eigenvalue equation 



M0„ = Xn(Pn (26) 

was solved on a lattice where M is a (8 ■ L^) x (8 ■ L^) matrix. The volume dependence of the 
lowest eigenvalues was studied and it was found that c ~ 0.5. To obtain a and b we notice that 
the quantities A„ and F^/ A^ were obtained for the lowest 292 nontrivial eigenvalues. In Fig. (1) 
we display both n and Fn/Xl^ vs. A„ on a log-log plot for L = Since n oc dX' p{X'), 
the slopes of these two curves can be used to determine a and b. In Fig. (2) we show the 
residuals (actual values divided by fit values) of the linear fits, where the least squares weights 
are determined by the statistical errors]^ We also perform fits for the L = 12 and 16 data. 
The structure observed in these residuals, the oscillations and the increasing departures for 
the smallest A„'s, are due to finite volume effects. The observed oscillations can be expected 

^An earlier lattice study lOJ for SU{2) found a sa .25 and b « .62. 

"^The authors of [11 kindly provided me with their raw data for the L — 12, 16 and 24 cases. 
^Very similar results follow from an unweighted fit where the lowest 27 eigenvalues are dropped. 
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to continue for some range of n > 292 as well. This makes the extracted slopes somewhat 
sensitive to the range of n used, and so these results for n < 292 should only be considered as 
indicative. 




1 b . ^ ^ 

0.01 0.1 

K 

Figure 1: The quantities n and Fn/\\ are plotted vs. A„ for the lowest 292 eigenvalues for 
L = 24, from a lattice study pT] . 



1.1 





(FXVfit 
n/fit 



















0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 

K 

Figure 2: The residuals of the linear fits (actual divided by fit values). 

We display our fits for a or 6 in Table 1, and the d values assume c = 0.5. We see that 
the L = 12 and 16 results do not satisfy the constraint a + b > 1/2. Although the results 
in Table 1 have uncertainties due to the finite volume effects, the trend of decreasing d for 
increasing L is clear. Presumably larger lattice sizes are needed to find the limiting value of d. 
In the cosmological context a value of d close to unity may give a vacuum energy of the correct 
magnitude. The results in Table 1 are still consistent with this, and they clearly motivate 
more refined lattice determinations of the scaling exponents. It is interesting that cosmology 
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can provide at least a lower bound 
simulations. 



on a scaling 



parameter measured in lattice gauge theory 



L 


a 


h 


a + h 


d 


12 


.094 


.238 


.33 


2.34 


16 


.146 


.285 


.43 


1.96 


24 


.205 


.381 


.59 


1.33 



Table 1: Fit values for the scahng exponents appearing in (15), (23) and (25). 



3 The cosmology 

In the last section we obtained the shift in the QCD vacuum energy density Ap^ac due to some 
change of physics occurring beyond a length scale L. This introduced a dependence of Ap^ac 
on 1/L, and now in the cosmological context we consider the replacement of 1/L by H{t). 
Hubble expansion changes the physics on and beyond Hubble scales in the sense that the 
vacuum fields are not able to attain the true Minkowskian vacuum state of QCD. It appears 
reasonable that the Hubble scale is the relevant scale, since points separated by more than 
the Hubble distance are receding from each other faster than the speed of light. This is true 
whether or not the universe is accelerating. Certainly the Minkowski description of the QCD 
vacuum would break down on these scales, since the interactions and correlations on larger 
scales would be disrupted. In particular it is fair to assume that the enhancement that these 
long wavelength modes enjoy in the true QCD vacuum no longer occurs. This is all that is 



needed to support our discussion below (22) 



So assuming that this residual energy density exists, the question is how it affects the 
evolution of the universe. This energy density cannot be described by some local term in an 
action, constructed to reproduce the H{t) dependence, since such a term would respond to all 
wavelengths. Implementing an effect only on wavelengths > H{t)~^ through the action would 
be highly nonlocal. Nevertheless there should be an effective energy- momentum tensor T^, 
depending on H{t), that can serve to describe the back-reaction effects on the cosmological 
evolution. There is also no reason for such a to be proportional to g^y, since over the 
relevant distance scales the Lorentz violating nature of the preferred cosmological frame is fully 
manifest. Nevertheless, in the cosmological rest frame, should preserve spatial isotropy 
and thus it should have a perfect fluid form in terms of the energy density p d and an effective 
pressure pr,. 

Next is the question of whether the dark energy exchanges energy with the normal matter 
and radiation content of the universe. Since the typical wavelengths of these normal compo- 
nents are vastly smaller than the dark energy wavelengths, we do not expect energy exchange. 
If there are fields with modes of Hubble wavelengths being excited for other reasons, then we 
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can simply redefine to include these contributions (assuming that they give subdominant 
contributions). We shall thus suppose that this full dark energy does not interact with 
normal matter and radiation. In this case we must have conservation, d^Tj^^^ = 0. This will 
relate the effective pd to po, with both quantities vanishing with H. 

Let us consider the evolution of quantities in terms of the scale parameter a rather than 



time. The dark energy density is Pd{(i) = CH{aY with our estimate from (24) being C ~ 
AQ'^'^3/327r^. The Friedmann equation becomes 

p{a) + CH{aY = -^H{a)\ 0<d<2. (27) 

In the early universe the normal energy density p{a) starts off larger than Pd{cI'), but it 
decreases faster as a — )■ oo. In this limit H{a) is driven towards a constant 

H(a^oo)^('^f' (28) 

from above. Thus the universe is driven to a de Sitter expansion. 

When we insert this limiting value back into Pd(o) = CH{aY and require that this corre- 
sponds to the observed value of dark energy (2.4 x 10~^ eV)*^ (and correct for the fact that the 
latter is the present value rather than the limiting value), we obtain an estimate for d. If we 
take Aqcd = 300 MeV then d ~ 1.026. Meanwhile d = 1 and 1.05 correspond to Aqcd = 130 
and 650 MeV respectively. And one should keep in mind uncertainties in the estimate of C. 

It is useful to keep d in the formulas (some similar results for d = 1 have appeared recently 
in [H]). We can assume that p(a) describes only pressureless matter, so p{a) = po/a^, and we 



take a = 1 at the present time. From (27) we can obtain 



_ H{a) 30^(a) 

a d-2-dnUay ^ ' 

where f2m(a) is the matter fraction of the total energy density. Using this we can determine 
the inflection point where d{t) = 0, or in other words when {aH{a))' = 0. It occurs when 

2 - d 

^m(ainfl) = 1- (30) 

S — d 

Another dimensionless quantity of interest at the present time is to-f^o- can be obtained 
through a numerical integration of d{t) = a{t)H{a(t)). 

Values of the inflection point red shift, Zmfl = l/amn — 1, and to-f^o can be obtained as a 
function of the present value of Qm = ^m(l)- In Table 2 we compare the d = 1.05 and 1 values 
with the d = values of a cosmological constant cosmology (note the different Qm values). 
We see that ZinA and toHo for the two cosmologies become more similar if Qm is reduced by 
about 0.04 or 0.05 when going from = to ci = 1 or 1.05. 
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.838 


.998 


-.819 


.21 


.793 


.986 


-.812 


.22 


.749 


.975 


-.804 


.23 


.708 


.964 


-.797 


.24 


.669 


.954 


-.790 



d= 1.05 





^infi 






.20 


.857 


1.004 


-.833 


.21 


.811 


.992 


-.826 


.22 


.768 


.980 


-.820 


.23 


.728 


.970 


-.813 


.24 


.688 


.959 


-.806 



d=l 
Table 2 
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.24 


.850 


1.025 


-1 


.25 


.817 


1.014 


-1 


.26 


.786 


1.003 


-1 


.27 


.756 


.993 


-1 


.28 


.726 


.983 


-1 



d = 




Figure 3: The effective equation of state parameter for dark energy when Qm = 0.22. a = 1 
is the present. 

The conservation equation p'^ + 3(1 + WD)pD/ci = and (29) determine the effective 
equation of state parameter, 

, , daH'ia] 



WD[a 



3 H{a) 
d-2 



(31) 



2 + rf(^]„(a) - 1) 

Thus WD{a) evolves from a value of (c/ — 2)/2 when dark energy is negligible to a value of — 1 
when it comes to dominate. At the inflection point WD(a.infl) = (c? — 3)/3. The present values 
Wd = wd(1) are given in the Table 2 while examples of the evolution of ^^(a) are shown in 
Fig. (3)j^ And flnally the form of the basic equation that describes structure formation in 

^In the case d = 1 wc have the more explicit result 

\ -2 

a4 (1 - n„,f , a2 (1 - n„,) 



f2m(a)/4a 



/4a- 
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normal matter 6 + 2H6 — ^QmH^^ = does not change, where again the dark energy does 
not participate at sub-Hubble wavelengths. 

For the cosmology then, the main observation is that the present values Qm and wd 
are both reduced as compared to the cosmological constant cosmology. These reductions 
appear to be consistent with present data. In particular we highlight a recent analysis 
[15] of a combination of datasets that have a similar intermediate redshift, with the goal 
of reducing systematic uncertainties. The values as quoted there are flm = 0.22tQ;Qg and 
wd = — 0.8llo;}8(stat)±0.15(sys). The analysis of more recent and future datasets, as is al- 
ready underway, will presumably go further at distinguishing a QCD generated dark energy 
from a standard cosmological constant. 

An interesting implication of these considerations is that they constrain the existence of 
confining gauge theories with mass scales higher than Aqcd- Their contribution to dark energy 
could be too large unless their corresponding d parameter is also larger. Conversely if such a 
contribution was the dark energy, then the larger value of d would give a cosmology further 
from the cosmological constant cosmology. For example if A = 1 TeV then d ~ 1.245, which 
gives Wd = —0.73 for Q^n = 0.22. 

To sum up, QCD is confining and has a mass gap. Nevertheless the long wavelength 
structure of the QCD vacuum differs markedly from that of free fields, either massive or 
massless, as can be seen in its response to the introduction of a color charge. We have argued 
that this implies significantly enhanced long wavelength contributions to the vacuum energy as 
well. But the true Minkowskian QCD vacuum presumably cannot be attained in an expanding 
universe, and this leads to a residual energy density that can be discussed independently of 
the cosmological constant problem. The true vacuum is approached as the universe expands, 
but eventually the residual energy dominates and the universe becomes stuck in a de Sitter 
expansion. The Minkowski vacuum is never attained. The question of whether the QCD 
residual energy is of the right size to be the present dark energy remains to be seen. It is 
interesting that cosmological observations and lattice studies can both come to bear on this 
question. 
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